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ABSTRACT 

Recent studies have focused on the the role of initially weak toroidal magnetic fields embedded 
in a stellar wind as the agent for collimation in planetary nebulae. In these models the wind is 
assumed to be permeated by a helical magnetic field in which the poloidal component falls off 
faster than the toroidal component. The collimation only occurs after the wind is shocked at 
large distances from the stellar source. In this paper we re-examine assumptions built into this 
"Magnetized Wind Blown Bubble" (MWBB) model. We show that a self-consistent study of the 
model leads to a large parameter regime where the wind is self-collimated before the shock wave 
is encountered. We also explore the relation between winds in the MWBB model and those which 
are produced via magneto-centrifugal processes. We conclude that a more detailed examination 
of the role of self-collimation is needed in the context of PNe studies. 

Subject headings: ISM: jets and outflows — magnetic fields — magnetohydrodynamics: MHD 



1. INTRODUCTION 

For over 30 years now, a purely hydrodynamic 
theory of planetary nebulae (PNe) shaping, known 
as the "interacting stellar winds" (ISW) model 
has been studied in considerable detail. The basic 
idea behind the ISW model is a fast, tenuous wind 
which overtakes a slow, dense wind from an earlier 
red giant or supergiant evolutionary phase. Appli- 
cations of this interacting stellar winds model to 
elliptical or bipolar PNe relies upon a generaliza- 
tion (GISW) which postulates the existence of an 
equatorial density enhancement in the slow wind 
(Kahn & West 1985; Balick 1987; Icke 1988). 
This density enhancement could be attributed to 
a binary system at the origin of the out-flowing 
winds (Soker 1995; Rasio & Livio 1994). This 
is especially true if the binary system suffers com- 
mon envelope evolution. For isolated star systems 
other possible explanations for the equatorial den- 
sity enhancement exist. One physical phenomenon 
possible of producing slow winds with aspherical 
densities is described by the "Wind Compressed 
Disk" model of Bjorkman & Cassinelli (1993) in 
which streamlines from a sufficiently rapidly rotat- 
ing star converge on the equatorial plane. Another 



possible mechanism for producing enhanced equa- 
torial density distributions via an initially dipolar 
field is described in the recent paper by Matt et 
al. (2000). 

As successful as this GISW model has been for 
describing PNe, there exist some prominent fea- 
tures which appear in new, high resolution images 
that can not embraced with the purely hydrody- 
namic model. These features include the presence 
of point-symmetric knots or ansae and jets as well 
as complex multi-lobed outflow structures (Balick 
2000; Manchado et al. 2000; Sahai 2000). While 
the GISW model can produce narrow jets (Icke 
et al. 1992; Frank & Mellema 1996; Mellema & 
Frank 1997) a large scale gaseous torus may be re- 
quired. Thus the new features seen in observations 
make alternative models appear more attractive, 
especially for those nebulae with coUimatcd point 
symmetric structures. One of the most attrac- 
tive alternative models which has been proposed is 
based upon another generalization of the interact- 
ing stellar winds scenario as originally described 
by Chevalier & Luo (1994). In this process, 
which we will refer to as the "Magnetized Wind 
Blown Bubble" (MWBB) model, the fast wind 
carries a weak magnetic field which is strength- 
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encd after passing through the wind shock. The 
nebula is thereby shaped by the latitudinal varia- 
tion in the total pressure in the hot bubble. This 
model has been studied numerically by Rozycka 
& Franco (1996) and has been subsequently gen- 
eralized (Garcia-Segura 1997; Garcia-Segura et 
al. 1999) by including a latitudinal variation in 
the density and velocity of the slow and fast wind 
based upon the "Wind Compressed Disk" model 
of Bjorkman & Cassinelli (1993). 

In this paper we critically examine assumptions 
built into the MWBB model. We have studied the 
model both analytically and muncrically and we 
show that there is a restricted regime for which 
the assumptions built into the MWBB model can 
hold. This regime is what must be called the very 
weak field limit and we show that to date numeri- 
cal simulations which show collimation are outside 
of this limit. In §2 we give a brief description of 
the MWBB model for the sake of completeness 
and clarity. In §3 we describe some of the previ- 
ous numerical studies of the MWBB model and 
the results. In §4 we study the fast wind region 
appropriate for the MWBB model and show that 
outside of the very weak field regime the fast wind 
will be magnetically coUimated. The conclusions 
of this section are applicable to the MWBB model 
as well as its generalization (Garcia-Segura 1997; 
Garcia-Segura et al. 1999). In §5 we discuss the re- 
lationship between MHD wind acceleration, colli- 
mation as well as its relevance to the jets observed 
in proto-planetary nebulae (PPNe). We summa- 
rize the results of this paper in §6. 

2. THE MWBB MODEL 

The purpose of this section is to provide an 
overview of key ingredients in the MWBB model 
as originally presented by Chevalier & Luo (1994). 
The model is based upon the interacting winds pic- 
ture for planetary nebulae (for a review see Frank 
(1999)). It is assumed that the flow can be de- 
scribed by four regions separated by three sur- 
faces: the wind shock, the contact discontinuity, 
and the ambient shock in order of increasing ra- 
dius. 

Similar to the Parker (1958) solar wind solu- 
tion, the MWBB model posits that the unshocked, 

fast wind can be described as a spherically expand- 
ing flow with foot-points tied to a rotating stellar 



surface. Applying flux conservation to this quasi- 
hydrodynamic model results in a poloidal mag- 
netic field Bp which scales as 1 /r^ and a toroidal 
magnetic field which for r ^ R,, scales as 
1/r, where i?* is the stellar radius. For a radial 
wind with constant velocity Vyj the density scales 
as l/r^; hence, the ratio a = B"^/ {Anpv^) is con- 
stant in the freely expanding wind. This ratio, a, 
is the main parameter for the flow. 

Next consider the region between the wind 
shock and the contact discontinuity, the so-called 
hot bubble. It is argued that the shock front in the 
stellar wind should occur at ^ 10^^ cm. Using the 
scaling of the magnetic field with radius and rea- 
sonable values for the stellar rotation rate and fast 
wind velocity Chevalier & Luo (1994) argue that 
at radii of interest B^/Bp > 10^. Hence they ne- 
glect the dynamical influence of the poloidal com- 
ponent of the magnetic field in the hot bubble. It 
is further assumed that the total differential dv/dt 
can be neglected. Hence, the flow in the hot bub- 
ble is steady, i.e. locally independent of time, with 
streamlines given by straight lines in the poloidal 
plane. The benefit of this assumption is that the 
calculation of the gas pressure and toroidal field 
is decoupled from the velocity and the solution is 
readily obtained. 

The region between the contact discontinuity 
and the ambient shock is argued to be strongly 
radiating. Thus it can be described in the "thin 
shell approximation" . The slow wind was taken to 
be spherically symmetric since Chevalier & Luo 
(1994) wanted solely to study the effect of the 
anisotropic (magnetic) pressure in the hot bubble 
in shaping the nebula. 

Using a formulation by Giuliani (1982), Cheva- 
lier & Luo (1994) derived and analyzed a set 
of coupled time dependent equations for the sys- 
tem. By assuming self-similarity they were able 
to remove the time dependence from the equa- 
tions. The resulting set of equations still requires 
numerical integration, though a limiting case for 
small parameters is analyzed analytically. The 
result of this limiting case suggests that devia- 
tions from spherical symmetry for the termina- 
tion shock can be expected for a > 0.006. This 
is further supported by the results of their numer- 
ical calculations. To put numbers to this, consider 
the situation for the PPN phase with a = 0.006, 
M = 10~^ Mq/uv, and Vu, = 300 km/s at a radius 
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of r 



10 



17 



This gives a toroidal magnetic 



field strength of i?^ 10 fiG. (Note that in this 
example we have used slower wind speeds which 
are appropriate to the PPN phase.) 

3. NUMERICAL STUDIES 

The first numerical simulations of the MWBB 
model were carried out by Rozycka & Franco 
(1996). Following Chevalier & Luo (1994) they 
neglected the poloidal component of the magnetic 
field and included only a toroidal magnetic field in 
the fast wind. Their simulations were cylindrically 
symmetric, on a grid with 300 x 500 cells and a 
resolution of about 10^'^ pc in each direction. In 
their study they found that the resulting flow was 
fairly insensitive to the angular variation in the 
toroidal magnetic field strength. They also found 
that for CT < 0.05 the flows evolve in a homologous 
manner as one would expect for a purely hydro- 
dynamic ISW model. Only for a > 0.05 did the 
simulations develop a strong pole- ward flow in the 
hot bubble. The authors ascribe this phenomena 
to the magnetic tension associated with the shock 
strengthened toroidal magnetic field. The numer- 
ical and analytic studies of the MWBB model dis- 
agree on the minimum value of a necessary for 
magnetic shaping of the hot bubble. Rozycka & 
Franco (1996) suggest that for a < 0.05, the 
length scale may simply be much larger than the 
simulation domain. Other studies of the MWBB 
model (Garcia-Segura 1997; Garcia-Segura et al. 
1999) included the presence of a gaseous torus as 
in the GISW model. In these studies a limiting 
value of a > .01 was found for effective magnetic 
collimation. 

We note that in the MWBB model the freely 
expanding wind has an unbalanced Lorentz force 
associated with the helical magnetic field. This 
force has a component perpendicular to the wind 
velocity. Hence, on physical grounds, it is reason- 
able to expect that on some length scale and for 
some range of field strengths, or cr, the fast wind 
will coUimate on its own. The fast wind may ex- 
perience significant collimation before it interacts 
with the wind shock. The range of initial parame- 
ters over which this occurs sets the limits on which 
the fast wind collimation can be reasonably ne- 
glected. It is to this question which we next direct 
our attention. 



4. FAST WIND COLLIMATION 

A basic premise of the magnetized interacting 
stellar wind model of Chevalier & Luo is that the 
fast wind velocity will remain spherically symmet- 
ric and that the will scale as 1/r for large r. 
This assumption can also be seen in the numerical 
simulations by noting that the radius of the "wind- 
sphere" on which the fast wind boundary condi- 
tions are set is ~ 10^^ — > 10^^ cm. The Lorentz 
force, however, is unbalanced, with components 
in a direction orthogonal to the velocity vector. 
While weak fields can not appreciably change the 
magnitude of the momentum, they can deflect the 
flow over large length scales. The question which 
remains to be answered is, "will the fast wind col- 
limate on length scales of interest and how does 
this collimation depend upon the parameter cr?" 

We approach this problem from a perturba- 
tion analysis in spherical polar coordinates (r, 0, </>) 
with the stellar rotation axis aligned with the ra- 
dial direction for = 0. We will present two cal- 
culations of the deflection, both of which assume 
purely radial flow for the zeroth order solution. 
The first calculation is presented as a simple, con- 
crete example and is based upon the Parker wind 
as described by Brandt (1970). This wind model 
treats the wind as purely hydrodynamic and ne- 
glects the dynamical influence of magnetic fields 
and rotational motion. The second calculation ad- 
dresses concerns about the appropriateness of the 
analysis for small r and moderate to large field 
strengths and is based upon the (Weber & Davis 
1967) model as described by (Barnes 1974). As we 
will see, both calculations are in good agreement. 

4.1. THE PARKER WIND 

The Parker stellar wind model is characterized 
by the following set of assumptions. The gas pres- 
sure, density and wind velocity are spherically 
symmetric {ve = = 0). The magnetic field 
is given by 

Br = -B* ( 

r 



Be = 
Bs = —Br 



rOsin(6') 



r 



where is the stellar radius (the base of the stel- 
lar corona), f2 the angular velocity, and the 
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poloidal component of the magnetic field at the 
stellar surface. For simplicity we take to be a 
split monopole field and assume solid body rota- 
tion. 

To analyze the collimation we need to solve 
Newton's equation for the motion in the 6- 
direction. For Vr ^ vg, Newton's equation of 
motion in the ^-direction is approximated by 



pVr d 
r or 



ee--{3 X B) 

c 



For ^ 7r/2 the Lorentz force in the ^-direction is 
given by 

, -B^cot^ 



2'Kr 

Since wc seek a closed form solution for the 
wind collimation, rather than study any partic- 
ular model for the wind acceleration we note that 
for r > i?« 



\Ba,\ > Br 



1 - 



where Voo is the asymptotic value of the radial 
velocity and we take Br and il as positive. Us- 
ing this expression for the toroidal magnetic field 
strength we will obtain a lower bound on the de- 
flection. The equation of motion in the ^-direction 
thus becomes 



-M = -2sm^cos^(^^^^^ 



r 



Next we note that the Michel velocity Vm (Michel 
1969; Belcher & MacGregor 1976) given by 



ATTpVr 



is a constant for steady, radial flow. Furthermore 
we note that following the definition of the param- 
eter a we have a = (Kn/Ko)^- Collecting terms 
we can write 



d / R 

— (rvg) = (-2aV^sm9cos9) ( 1 

or \ r 



The integration of this equation is straight for- 
ward. Applying the boundary condition that ve = 
at r = ii* we obtain 

v$ = —2aVoo sin ^ cos ^ (l — a;^ -|- 2xln(a;)) 



where x = R^^/r. Now to study the deflection or 
collimation of the fast wind we need to calculate 
the streamlines of the flow. That is, we wish to 
integrate the equation 



d^ 
dr 



V0 

rvr 



Once again wc content ourselves with a lower 
bound for the deflection and substitute Voo for Vr 
in this expression. This equation can then be in- 
tegrated, resulting in a relation which describes 
the path taken by a particle as it moves from the 
point (i?*,^*) on the stellar surface to the point 
(r, 6) some distance away. The solution is 



tan^ 
tan^* 



^2^+4ax exp{CT(5 -4x- x^)} 



where again x = R^/r. For r ^ -R* this equation 
is well approximated as 



tan 
tan 0* 



R, 



2a 




0.04 0.06 
a 

Fig. 1. This plot shows how the collimation of 
the fast wind depends on the parameter a for two 
radii, (r/i?*) = 10^, 10^ 

This calculation shows that the fast wind will 
collimate for sufficiently high values of a and long 
propagation lengths. In Figure 1 we plot the angle 
6 (in degrees) to which the gas is deflected as a 
function of a for 6*, = 45° and (r/i?*) = 10^ 10^. 
For example, for ii* = 10^^ cm Figure 1 shows 
that by the time gas parcels which originate at 
^* = 45° reach a radius r = 10^^ cm they will 
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be deflected to 6» w 35°, 22° for a = 0.02, 0.05 
respectively. 

We note that the Parker wind model is essen- 
tially a hydrodynamic wind model with the mag- 
netic field assumed weak enough to have no dy- 
namical influence. Thus one might worry that the 
preceding argument is only applicable for cr <C 1. 
Furthermore, it is well known that for radii less 
than or equal to the Alfven radius, Ra, the mag- 
netic field enforces near co-rotation of the plasma 
with the star. Thus there exits the concern that 
by assuming i;^ = we have overestimated the 
strength of the toroidal magnetic field at small 
radii and, thereby, overestimated the degree of 
collimation. In what follows we demonstrate that 
collimation also occurs when the perturbation cal- 
culation includes the plasma rotation (Weber & 
Davis 1967; Barnes 1974). 

4.2. THE WEBER-DAVIS WIND 

In this next calculation we will make use of the 
Weber & Davis (1967) formalism. As described by 
(Barnes 1974) wc assume! that vg = so that the 
results can be extended away from the equator. As 
in the previous calculation we will assume a split 
monopole field geometry, solid body rotation and 
spherically symmetric gas pressure and density for 
simplicity. Furthermore, we will assume that in- 
terior to the Alfven radius that vg = remains a 
good approximation. We will not assume that the 
rotational velocity, v^, is zero, but instead include 
it self consistently with the toroidal magnetic field. 

To analyze the collimation we solve Newton's 
equation for the motion in the 6'-direction. For 
Vr Vg, Newton's equation of motion in the 9- 
direction is approximated by 



r or 



pv\ cot 9 1 

=e,--(JxB) = F/'. 

r c 



For ^ 7^ 7r/2 the Lorentz force in the ^-direction is 
given by 

-Blcot9 

= 



27rr 

The equation of motion in the ^-direction can thus 
be rewritten as 

d , , - cot 9 I Bl \ 



Integrating the azimuthal component of the induc- 
tion equation gives 

B(i,Vr — Br{v^ — rfisin^) = 

where again f2 is the angular velocity at the base 
of the stellar corona. Integrating the azimuthal 
component of Newton's equation of motion gives 

BrB^ L{9) 



where L is the specific angular momentum. At the 
Alfven point (r = Ra) determined by the condi- 
tion Airpv'^ = B^., these two equations are linearly 
degenerate giving the so-called Alfven regularity 
condition, L = i?^f2sin^. Using these relations 
and a bit of algebra we can solve for and B^ 
in terms of parameters evaluated at the Alfven ra- 
dius. The equation of motion in the ^-direction 
can then be rewritten in the form: 



d_ 

dr 



{rvg) = 



-sine cos f2V: 



:(l_y)^_r2n2y2 



where Vm is the Michel velocity (defined in the last 
section) and 



Y = 



{RA/rf - pIpa 

1 - pIpa 



The parameter y is a measure of the degree of 
co-rotation of the gas due to the torque applied 
by the magnetic field and rotational motion of the 
magnetic foot-points. That is, F = describes 
the situation of no rotational motion, = 0, 
and Y = 1 describes the situation of perfect co- 
rotation, Bfj, = 0. The value of Y evaluated at the 
Alfven radius, Ya, must be understood in terms of 
I'Hopital's rule and shows that Ya is also a mea- 
sure of the radial acceleration of the wind. 

For r > Ra, Y oc (Ra/t)^- Thus the equation 
of motion in the ^-direction shows that at large 
radii the force due to the magnetic field will dom- 
inate over the centrifugal force. Unfortunately, it 
is difficult to proceed further without some knowl- 
edge of Vr since this controls the radial variation of 
the parameter Y . One can show that for r > Ra, 



Y<Ya 



Ra 



when the radial velocity variation is slow enough 
that 



dvr 
dr 



< 



2{vr - Va)r: 



r{r-^ - R\) 
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It is under this condition that wc can calculate a 
lower bound on the coUimation of the wind with 
radius. The equation of motion in the ^-direction 
can be rewritten as: 



d_ 
dr 



(rve) 



-2crVoo sin cos { 1 — k 



R\ 



where 



Following the same procedure described for the 
previous calculation we obtain a relation for the 
streamline connecting the point {Ra,Oa) to the 
point (r, 9) which describes a lower bound on the 
collimation. The result is 



tan6' 
tan 9 a 



x^" exp{£7(2 + K - 2(1 + k)x + Kx^)} 



where x = Ra/t. For r » Ra this equation is 
well approximated as 



tan(? 
tan^A 



.(2+«) 



1(7 



For comparison with the previous calculation 
we would like to know, order of magnitude, a rea- 
sonable number for the exponent (2 + k)u. Using 
the definition of k we have 



(2 + K)a 



td2 o2y2 

2{1 + YA)a+^^^^ 



2^2 



To estimate the first term we need an approxima- 
tion for Ya- As described earlier, the parameter Y 
is a measure of the co-rotation of the gas. There- 
fore, Ya is a measure of the force balance in the 
^-direction at the Alfven radius. Balancing the 
0-component of the LoKuitz forcc^ and centrifugal 
force at the Alfven radius gives Ya = (2 — \/2). 
To estimate the second term in the equation for 
(2 -|- K)a, there are two relevant limits which will 
be discussed more fully in a later section. In the 
so-called slow rotator limit the second term can be 
arbitrarily small. In what is called the fast rota- 
tor limit, Vm « Voo and Va « (2/3) K„ (Belcher & 
MacGregor 1976). Evaluating the Michel veloc- 
ity at the Alfven Radius one finds = R\Q'^Va- 
Thus in the fast rotator limit we have 



(2 + k)(j ■ 



2(l + y^)a + Y 



Inserting Ya = (2 — \/2), wc find the exponentials 
in both this and the previous calculation to be 
of order ~ 1. Thus we find good agreement for 
the intrinsic collimation in both the Parker and 
Weber-Davis model calculations. 

Numerical simulations which start at large radii 
and assume a spherical flow are thereby con- 
strained to small values of cr. In order to study 
models with appreciable field strengths (fields be- 
yond the " very weak limit" ) requires following the 
flow evolution at small radii. This in turn neces- 
sitates including the poloidal magnetic field since 
at small radii the poloidal and toroidal fields are 
comparable in strength. In summary, we conclude 
that a full treatment of the propagation and mag- 
netic collimation effects is warranted for the study 
of magnetized outflows such as the MWBB model. 

5. WIND ACCELERATION AND COL- 
LIMATION 

The collimation of magnetized winds has been 
extensively studied. In most investigations, how- 
ever, both the collimation and acceleration (also 
know as launching) are considered together (Pu- 
dritz & Konigl 2000). In view of the rather ex- 
tensive literature on the subject it is useful to at- 
tempt to put the MWBB model into the context 
of these investigations. This is particularly true 
because jets often occur in the PPNe phase (Sa- 
hai & Trauger 1998; Sahai 2000) before the star 
has become hot enough to radiatively drive a fast 
wind (Frank 2000). The launching of the wind is 
not addressed in the MWBB model. 

The work of Weber & Davis (1967) was the 
flrst to explore the role of the magnetic field in 
launching a wind. This work was carried out in 
1.5-D, focusing only on the launching of a wind 
at the equator of a rotating star. The ability of a 
magnetized rotator to launch and collimatc a wind 
was explored in axisymmetric calculations in the 
work of Blandford & Payne (1982) and Sakurai 
(1985). Since that time many workers have ex- 
plored magneto-centrifugal winds (see Lamers & 
CassineUi (1999) or Pudritz & Konigl (2000) for 
a recent review of the subject). We note that a 
majority of the work on wind launching and colli- 
mation has focused on the role of accretion disks 
as the source of the wind. Recently, however, 
their has been some renewed interest in a star as 
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the magnetized rotator source of the the magneti- 
cally launched and collimated wind (Bogovalov & 
Tsinganos 1999; Tsinganos & Bogovalov 2000). 

A fundamental assumption of the MWBB 
model is that the magnetic field in the wind is 
sufSciently weak that intrinsic MHD wind coUima- 
tion can be neglected. We now address the issue 
of the self-consistency of this assumption in terms 
of the acceleration of the wind. Belcher & Mac- 
Gregor (1976) utilized the Weber-Davis model to 
study the transition from stellar winds in which 
the magnetic field is insignificant to those in which 
the magnetic field is dominant in driving the wind. 
The former and later were classified as slow and 
fast magnetic rotators. This classification can also 
be described by comparing the Michel velocity Vm 
to the Parker velocity Vp, where Vp is defined as 
the velocity at infinity in the absence of rotation. 
Slow magnetic rotators have V„i ^ Vp, while fast 
magnetic rotators have Vm ^ Vp. The important 
point to note here is that fast magnetic rotators 
drive winds via the transfer of Poynting flux into 
kinetic energy flux. In the fast magnetic rotator 
limit Voo ^Vm- 

Consideration of the parameters used in MWBB 
models shows that the majority of the MWBB 
simulations carried out to date (with a < 0.1) 
are not in the fast rotator regime. (Recall that 
c = {Vm/Voo)^ ■) Thus these models arc consistent 
with the assumption that the wind is launched by 
some means other than conversion of Poynting 
flux to kinetic energy flux. Conversely the winds 
used in MWBB studies can not be said to have 
come from a form of magnetized disk/star rotator 
wind model since these should have Vm/Voo ~ 1- 
This is an important point as it speaks to the ori- 
gin of the wind in the PPNe stage. The results of 
our perturbation analysis are of interest because 
they show that even a field that can be considered 
weak {i.e. a <C 1) can produce a pre-coUimation 
in the wind before the wind is shocked and the 
MWBB mechanism becomes operative. 

Finally we note that other studies of magnetic 
winds from slow rotators have shown some degree 
of coUimation. Tsinganos & Bogovalov (2000) 
found that the solar wind could achieve a signifi- 
cant degree of collimation at large distances when 
the solar rotation rate was increased by only a 
factor of 5. These results are consistent with our 
perturbation analysis. 



6. CONCLUSIONS 

In this paper we've presented an analysis of 
the MWBB model. We've shown that outside 
of the "very weak field" limit, the inherent mag- 
netic collimation of the fast wind is neglected. In 
this sense, the numerical models which begin at 
large radii and include magnetic fields beyond the 
very weak field limit are not consistent with the 
assumed conditions at small radii. In addition 
we have explored the issue of intrinsic collima- 
tion in the MWBB winds in light of the combined 
launching and collimation of winds that is possible 
through magneto-centrifugal processes (Tsinganos 
& Bogovalov 2000). An analysis of conditions 
used in MWBB studies shows that these winds 
can not originate from stars lying in the fast mag- 
netic rotator regime. This point is relevant in de- 
termining the applicability of the MWBB model 
PPNe. Many PPNe show collimated flows even 
though their central stars do not have the photon 
momentum flux required to drive these outflows 
(Alcolea et al. 2000; Sahai & Trauger 1998). The 
role of magnetic rotators in driving PPNe flows 
has recently been articulated by Blackman, Frank 
& Welch (2000) and Blackman et al. (2001). 

We wish to emphasize, however, that the funda- 
mental processes described by the MWBB model 
are sound and, in fact, represent a class of behav- 
ior that has not been given adequate attention in 
the literature on MHD collimation. As the pre- 
vious MWBB studies have shown the post-shock 
hoop stresses in a magnetized wind with signifi- 
cant toroidal field will act as an additional colli- 
mation mechanism above and beyond the intrinsic 
collimation provided via the launching process. 

Thus magnetohydrodynamic collimation and 
shaping of Planetary nebulae as well as produc- 
tion of such objects as the knots and jets observed 
within PNe is a very viable mechanism. The great 
promise which lies in the incorporation of mag- 
netic fields into models of PNe formation is in the 
description of the shaping mechanisms. With this 
should also come a better understanding of the 
period of stellar evolution during which the mass 
loss is evolving from a slow to fast stellar wind. 
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